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ABSTRACT 

Locating the centers of dark matter halos is critical for understanding the mass profiles of halos as 
well as the formation and evolution of the massive galaxies that they host. The task is observationally 
challenging because we cannot observe halos directly, and tracers such as bright galaxies or X-ray 
emission from hot plasma are imperfect. In this paper we quantify the consequences of miscentering 
on the weak lensing signal from a sample of 129 X-ray selected galaxy groups in the COSMOS field 
with redshifts < z < 1 and halo masses in the range 10 13 — 10 14 M Q . By measuring the stacked 
lensing signal around eight different candidate centers (such as the brightest member galaxy, the 
mean position of all member galaxies, or the X-ray centroid), we determine which candidates best 
trace the center of mass in halos. In this sample of groups, we find that massive galaxies near the 
X-ray centroids trace the center of mass to < 75 kpc, while the X-ray position and centroids based on 
the mean position of member galaxies have larger offsets primarily due to the statistical uncertainties 
in their positions (typically ~ 50 — 150 kpc). Approximately 30% of groups in our sample have 
ambiguous centers with multiple bright or massive galaxies, and some of these groups show disturbed 
mass profiles that are not well fit by standard models, suggesting that they are merging systems. We 
find that halo mass estimates from stacked weak lensing can be biased low by 5 — 30% if inaccurate 
centers are used and the issue of miscentering is not addressed. 

Subject headings: cosmology: observations - galaxies: clusters - gravitational lensing: weak 



1. INTRODUCTION 

Galaxy groups and clusters are important sites of 
galaxy evolution and the abundance of these massive 
objects provides a sensitive probe of the amplitude of 
matter fluctuations and other cosmological parameters. 
Analyses of these structures require some knowledge of 
the location of the centers of their gravitational poten- 
tials. Because the total mass distribution is dominated 
by dark matter and is not directly observable, halo cen- 
ters are typically assumed to be traced by a massive 
galaxy or the density peak of radiating hot gas. Mis- 
centering is a critical issue when estimating the masses 



of groups and clusters because it adds significant system- 
atic uncertainties (e.g., Johnston~et al. | 2007a|b Mandel- 
baum et al. 2010 Rozo et al. 201ip , and also de grades 



constraints on the co ncentration of mass profiles (Man 



delbaum et al. 2008). Velocity offsets between observa- 
tional tracers and halo centers impact studies of satellite 



kinematics (Skibba et al. 2011 Wojtak et al 



redshift-space distortions (Hikage et al 



2011) and 



2012). On the 
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other hand, offsets between observational tracers and the 
true halo centers can provide information about the dy- 
namical state o f these systems and about the properti es 
of dark matter ( |Clowe et aL]|2006| |Massey et aT]|2011| ). 

Finding halo centers is challenging for a number of rea- 
sons. Galaxy formation models (as well as halo models 
for describing the multiplicity of galaxies within halos) 
typically place the brightest or most massive galaxy at 
the center of each halo. But the brighte st galaxy in a 
cluste r is not always the central galaxy ( |Skibba et al. 
2011 and references therein). Groups and clusters form 



from mergers of halos where the most massive halo be- 
comes the host halo with its central galaxy, and smaller 
halos become subhalos with satellite galaxies. Several 
analyses of data from group catalogs and field surveys 
have found that there is some intrinsic scatter in stel- 



lar mass and luminosity at fixed halo mass (Yang et al. 
20091 [More et~aLl E009l I Leauthaud elTaTpol^ |Reddick 
et al. 2012p , which implies that halos can end up with 
satellites that are intrinsically more massive or luminous 
than the central galaxy. Additionally, there are uncer- 
tainties in measuring any observable quantity such as 
stellar mass that can cause a satellite to be misidenti- 
fied as the most massive central galaxy, and structure 
projected along the line of sight can confuse the identi- 
fication of member galaxies. Another difficulty is that 
merging systems are dynamically unrelaxed, which can 
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produce offsets between the central galaxy and the halo 
center or other tracers such as the X-ray center. The 
systematics introduced by picking centers that do not 
coincide with the "true" center of mass are important 
and need to be quantified. 

Gravitational lensing is a powerful tool for finding the 
centers of mass of halos since it is sensitive to the total 
mass distribution along a line of sight. Mass maps can be 
constructed for individual systems with stro ng lensing or 



for massive clusters with weak lensing (e. g. , |Smith et al. 
2005] |Oguri eTaT1[20T0{ |Shan et al.|[2010| . Such studies 



otten rind reasonable agreement between the positions 
of bright massive galaxies, X-ray emission, and lensing 
mass peaks, with a handful of interesting examples that 
illustrate how dark matter peaks can be offset from hot 



gas in merging gal axy clusters (e.g., Clowe et al. 2004 
Bradac et aL|2008[ ) 



In this paper, however, we are concerned with a large 
statistical sample of groups with lower masses and higher 
redshifts, a regime of interest for many current and fu- 
ture weak lensing surveys. The typical signal-to-noise 
ratio for the weak lensing signal from individual groups 
is low, so we cannot identify their halo centers individ- 
ually. By stacking the lensing signal from many groups, 
we determine the mean mass profile around a given cen- 
ter. We repeat this process for different candidate centers 
and compare the resulting profiles to find the best tracer 
of the center of mass. The center of a smooth halo can 
be identified as the position where the lensing signal is 
maximized on small scales. Other components such as 
a massive galaxy or subhalo that is offset from the halo 
center could produce an additional peak in the lensing 
signal, so we must account for that when modeling the 
signal. 

We analyze a sample of 129 X-ray selected galaxy 
groups at redshifts < z < 1 fro m the COSMOS fiel d 
(fScoville et al.||2007| ), described in |George et~aT1 ( |2011[ ). 
With COSMOS data, we have X-ray luminosities and 
centroids for each group, with member galaxies identi- 
fied using photometric redshifts derived from over thirty 
ultraviolet, optical, and infrared bands, and a subsample 
with spectroscopic redshifts. We use weak lensing mea- 
surements from high resolution Hubble imaging to study 
the accuracy with which tracers such as bright galaxies 
and X-ray emission identify the centers of dark matter 
halos. 

This paper is the second in a series studying t he galaxy 



content of X-ray groups. |George et al. (2011 hereafter 
Paper I) presented a catalog of group membersnip assign- 
ments from a Bayesian treatment of photometric red- 
shifts, along with extensive tests of the selection algo- 
rithm using mock catalogs and subsamples with spectro- 
scopic redshifts. Initial analyses of group members were 
used in that paper to demonstrate an environmental in- 
fluence on galaxy colors out to z = 1. A previous weak 
lensing study of this group sample was used to constrain 
the m ean relation between X-ray luminosity and halo 
mass ( jLeauthaud et al.||2010 ). 

In this paper we study the centers of groups in de- 
tail to optimize observational choices for centering, to 
study the impact of miscentering on measurements of 
halo properties, and to explore the effects of merging 
and substructure on lensing measurements. The outline 
of the paper is as follows. Section [2] describes the data 



used in our analysis, including the X-ray group catalog, 
assignment of member galaxies, and lensing shape mea- 
surements. We define eight candidate group centers in 
Section |3j and describe our procedure for testing differ- 
ent choices of centers in Section |4j Section IS] presents 
the results of our analysis, and Sections [6] and]?] provide 
discussion and conclusions of our work. 

We adopt a WMAP5 ACDM cosmology wi th O m = 
0.258, Q A = 0.742, H = 72 h 72 km s^ 1 Mpc" 1 dDunkley 



2009 ) following the in i tial le nsing analysis of these 
groups by Leauthaud et al. ( |2010[ ). All distances are ex- 
pressed in physical units with h = 0.72. X-ray luminosi- 
ties are expressed in the 0.1-2.4 keV band, rest-frame. 
All magnitudes are given on the AB system. To approx- 
imate the virial radii of halos, we use i?200c which is the 
radius within which the mean mass density equals 200 
times the critical density of the Universe at the halo red- 
shift, p c (z). The corresponding mass enclosed within this 
radius is M 2 oo c = 200p c (47r/3)i?2ooc- We also assume ha- 



los follow a Navarro-Frenk- White (NFW, jNavarro et al. 
1996 1 density profile, with a concentration parameter 
C200C and a scale radius R s = -R2ooc/c2ooc- We use the 
term "group" to denote a set of galaxies occupying a 
common halo, and the halo masses of these groups is in 
the range 10 13 — 10 14 M n as estimated with weak lensing 
( Leauthaud et al.||2010 ). We will generally refer to more 
massive structures as clusters following convention, but 
make no other physical distinction between groups and 
clusters. 

2. DATA 

To study how the constituents of galaxy groups trace 
the centers of mass of dark matter halos, we use an X- 
ray s elected sample of ga laxy groups from the COSMOS 
field ( Scoville et al.|2007| ). We refer the reader to Paper I 
for details of the data and methods used to construct the 
group catalog as well as tests of its properties with simu- 
lations and spectroscopic data. In this section, we briefly 
describe aspects of the catalog that are relevant for cen- 
tering including the assignment of member galaxies to 
groups. We also introduce the galaxy shear catalog used 
in our weak lensing analysis. 

2.1. X-ray Group Catalog 

Our sample of galaxy groups has been selected from an 
X-ray mosaic combin ing images fro m the XMM-New ton 
( |Hasinger et al.|2007[ ) and Chandra ( [Elvis et al.|2009[) ob- 
servatories f ollowing the procedure of |Fmoguenov et al.| 
( |2009| [2010| ) . A wavelet filtering of theX -ray mosaic is 
used to distinguish extended structures on scales of 32" 
and 64" from contaminants on smaller scales like active 
galactic nuclei (AGN). Once extended X-ray sources are 
detected, a red sequence finder is employed on galaxies 
with a projected distance less than 0.5 Mpc from the cen- 
ters to identify an optical counterpart and determine the 
redshift of the group, which is then refined with spectro- 
scopic redshifts when available. 

A quality flag (hereafter XFLAG) has been assigned 
to each group based on the reliability of the opti- 
cal counterpart identification. We study groups with 
XFLAG=1 or 2, indicating a confident spectroscopic asso- 
ciation, while higher values indicate uncertain counter- 
parts which could be due to projection effects or photom- 
etry contaminated by bright foreground stars. Sources 
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with xflag=1 also have clear X-ray centroids, with an 
uncertainty in each position coordinate, ax, equal to the 
wavelet scale of 32" divided by the signal-to-noise of the 
flux measurement, while sources with xflag=2 have less 
certain X-ray centroids for which we have ax = 32" set 
by the wavelet scale of the flux measurement. 

To ensure a clean sample of groups with robust mem- 
bership assignment, we employ the additional quality 
cuts suggested in Paper I, excluding groups that are near 
field edges or have significantly masked areas, potentially 
merging groups identified as distinct X-ray sources but 
with significantly overlapping volumes, and poor groups 
with fewer than four members identified. After these 
quality cuts, we have 129 groups in our sample ranging 
from redshift < z < 1. 

2.2. Galaxy Membership 

To determine how well galaxies trace the matter distri- 
bution in groups, we must first identify the galaxies that 
reside in them. The COSMOS field has extensive imag- 
ing in over thirty u ltraviolet, optical, and infrared bands 



( Capak et aLp007 1 , enabling the determination of stellar 
masses (s ee Paper I for det ails) and precise photometric 
redshifts (jllbert et al.|2009| and Paper I for further tests) . 
In Paper 1, we presented a catalog of member galaxies for 
these X-ray groups, selected according to their photomet- 
ric redshifts and proximity to X-ray centroids. Briefly, 
a Bayesian membership probability, P mC m, is assigned 
to each galaxy by comparing the photometric redshift 
probability distribution function to the expected redshift 
distribution of group and field galaxies near each group. 
From the list of members with P mcm = 1 — Pfieid > 0.5, 
the galaxy with the highest stellar mass within an NFW 
scale radius of the X-ray centroid (including the posi- 
tional uncertainty, ax) is selected as the group center. 
We call this object the MMGG sca i c , for "most massive 
group galaxy within a scale radius" . A final membership 
probability is assigned by repeating the selection process 
within a new cylinder re-centered on this galaxy. 

We have extensively tested this selection algorithm us- 
ing mock catalogs and with subsamples of galaxies for 
which we have spectroscopic redshifts, and found it to 
be both pure and complete near group centers; within 
0.5i?2ooc and down to our limiting selection magnitude 
(F814W=24.2), 84% of galaxies selected as group mem- 
bers truly belong in groups, and 92% of true group mem- 
bers are selected as such. In this paper we use the mem- 
ber catalog derived from photometric redshifts, which 
has an average of 26 members per group. From that list, 
there are an average of 6 members per group with spec- 
troscopic redshifts for calibration and determining group 
redshifts. 

2.3. Weak Lensing Data 
The galaxy shape measurem ents used for our wea k 



lensing analysis are described in Leauthaud et al. (2007). 



These are derived from high-resolution imaging over 
1.64 degrees 2 of the COSMOS field with the Advanced 
Camer a for Surveys (ACS) on t he Hubble Space Telescope 
(HST; Koekemoer et al. 20071 to a limiting magnitude 
of F814W=26.4. Variations m the point-spread func- 
tio n (PSF) with positio n and time are modeled follow- 
ing Rhodes et al. (2007), and galaxy shapes are derived 



using the RRG method ( |Rhodes et aLpOOO) ). The PSF- 
corrected shapes are converted to estimators of shear, 7, 
using a shear susceptibility factor calculated from mo- 
ments of the global distribution of shapes and a calibra- 
tion factor determined from simulated images. Updates 
to the procedure a nd the shear catalog are described in 



detail elsewhere, in Leauthaud et al. (2012). These im- 



provements include a mor e detailedcorr e ction of charge 



transfer inefficiency from Massey et al. (2010), and an 
empirical derivation of the dispersion in shear measure- 
ments in bins of magnitude and detection significance. 
This estimate of the shear dispersion includes contribu- 
tions from intrinsic shape noise and shape measurement 
uncertainties, and varies from a 1 m 0.25 for bright galax- 
ies to <t 7 m 0.40 for faint objects. 

The stacked weak lensing signal is derived from the 
average tangential shear, j t (R), of background source 
galaxies at a projected distance R from the center of 
each group. The she ar is related to the ex cess surface 
mass density AE(i?) ( |Miralda-Escude|| 1 99 1 [ ) 

(1) 



AE(i?) = £(< R) - E(i?) = 7t (i?)S C rit, 



where E(< R) is the mean surface density within radius 
R and E(i?) is the azimuthally averaged surface density 
at R. The critical surface density E cr it is a function of 
the angular diameter distances between the observer (O), 
lens (L), and source (5), 



E r .rit — 



c 2 Dos 
4ttG DqlDls ' 



(2) 



where c is the speed of light and G is the gravitational 
constant. 

In order to compute E crit from Equation ([2]), we need 
to know the distances to both the sources and lenses. 
The group catalog provides lens redshifts which come 
primarily from sp ectroscopic data including zCOSMOS 
( |Lilly et ' al.|2009[ and in prep.). F or background sou rces, 
we use photometric redshifts from Ilbert et al. ( 2009[ ). To 
avoid contamination due to uncertainties in photometric 
redshifts, we use only sources with zs~~zl > max[0.1, a z ] 
where a z is the 68% uncertainty in the source redshift. 
We also exclude sources with a secondary p eak in their 
redsh ift density function (i.e. ZP2 7^ in the Ilbert et al. 
2009 catalog) which have a significant fraction of catas- 
trophic redshift errors. With these cuts, the source cata- 
log contains 210, 015 galaxies with well-measured shapes 
and redshifts, providing a source density of 36 galaxies 
per arcmin 2 . 

To obtain a significant measurement of AE, we must 
combine the signal from many lenses and background 
sources. The combined measurement is a weighted sum 
over pairs of lenses i and sources j , 



AE = 



J crit,ij 



(3) 



where the weight Wy = (E cr itcr 7j ij) -2 is the inverse vari- 
ance of the measurement. We measure AE in annular 
bins from 20 kpc to 1 Mpc. Covariance between mea- 
surements becomes an issue on larger scales where back- 
ground sources can be paired with multiple lenses, but 
this is not significant over the scales we measure. Uncer- 
tainties in AE are determined from the inverse square 
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root of the sum of the weights of lens-source pairs. 
3. DEFINING CANDIDATE CENTERS 

The "center" of a galaxy group requires some defini- 
tion. There is ambiguity in centering even when con- 
sidering simulated dark matter halos; the mass centroid, 
most bound particle, and density peak can all be differ- 
ent because of asphericity and substructure. The choice 
of group and cluster centers in observational data sets 
is further limited by the available measurements. Here 
we review a variety of definitions of group centers and 
their relative advantages. Our aim is to use weak lens- 
ing to determine which candidates most accurately trace 
(on average) the centers of mass of dark matter halos. 
We will consider a variety of candidate centers and begin 
by studying the level of agreement between them. Our 
choice of candidate centers is meant to explore the range 
of options available for multi-wavelength data sets while 
using a simple set of rules for identification; however it 
is not an exhaustive list of possible centers. 

It is useful to separate these definitions into two broad 
categories. We call the first set "galaxy candidates" since 
they are centered on a single galaxy, and the second set 
"centroid candidates" which are defined for a spatially 
extended quantity like the galaxy density field or X-ray 
emission and are in general not centered on an individ- 
ual galaxy. Some centering algorithms take a hybrid 
approach, using the proximity of neighborin g members 
to ultimate ly select a luminous galaxy (e.g., Robotham 
et al.pOll" ), but we do not test those methods here. 

When identifying centers based on the galaxy content 
of groups, we select from galaxies with membership prob- 
ability ? mcm > 0.5, as described in Section [272] and Pa- 
per I. Though the list of members is defined around one 
of the candidate centers (the MMGG sca i G ), the radius 
(-R200c) used to select members is large enough that the 
initial choice of center should not impact our results. 
Each of the centers based on galaxy fluxes (e.g., brightest 
group galaxy) use the observed magnitude in the F814W 
band, taken with the ACS on HST, with no corrections 
for dust or evolution. Since these measurements do not 
account for the change in rest-frame wavelength probed, 
they will be more sensitive to recent star formation at 
higher redshifts. Centers based on galaxy masses use the 
full measured spectral energy distribution (SED) so these 
effects are diminished. 

3.1. Galaxy Candidate Centers 

Many clusters have a central dominant galaxy with an 
extended stellar envelope, often located near the density 
peak of hot intracluster gas as seen in X-rays and the 
peak of the matter density pro bed by lensing or kine- 



or BGG i n groups), since the selection is relative ly sim- 



ple (e.g., Koester et al. 2007b Hao et al. 2010). But 



matics (e.g. 



and references therein). 



Lin fc Mohr||2004 
This motivates the choice ot a single galaxy to trace the 
centers of groups and clusters. The general picture is 
further supported by numerical simulations of dark mat- 
ter halos and subhalos, and has been encapsulated in the 
halo model which successfully describes many aspects of 
large-scale structure includi ng measurements of ga l axy 



je-g., 



Cooray fc Sheth||2002 



clustering and lensing 

havi et al. |2005| [Mandelb aum et al.||20Ubb| |LeauthaU' 
et al.||2012| ). 



1'hus a popular choice when defining cluster centers in 
optical catalogs is the Brightest Cluster Galaxy (BCG; 



the choice of filter and aperture used for the flux mea- 
surement has an impact on which galaxies are selected; 
differences in redshift, star formation history, and dust 
content affect the flux observed in a given band, so a 
single flux measurement cannot reflect the complicated 
physical processes occurring in group centers. Col or cuts 
can be used to is olate a few of these effects (e.g., Glad- 
ders fc Yee|2000[ ), though they often come with assump- 
tions about the properties of central galaxies. For exam- 
ple some group catalogs use the brightest red sequence 
galaxy to identify group centers, avoiding galaxies that 
are bright due to recent star formation in favor of massive 
galaxies with old stellar populations. 

Stellar masses are a promising alternative to observed 
or rest-frame luminosities since they correlate more di- 
rectly with the masses of halos in which galaxies reside. 
However, stellar mass estimates require more detailed 
measurements of the SEDs of galaxies and are fraught 
with larger uncertainties than simple fluxes or luminosi- 
ties. 

In this paper we consider four galaxy candidates, se- 
lected based on flux or stellar mass and distance to the 
X-ray position: 

• MMGGscaic) the galaxy within R s + ox of the X- 
ray centroid having the greatest stellar mass. 

• MMGGr200 : the galaxy having the greatest stellar 
mass of all group members within i?200c- 

• BGG sca i c : the brightest galaxy within R s + ox of 
the X-ray centroid. 

• BGGr2oo : the brightest galaxy of all group mem- 
bers within i?200c- 

The X-ray centroid (with uncertainty ctx) is used as 
the starting point for selecting these galaxies because it 
should roughly trace the halo center and we do not have 
lensing centers for individual groups. Note that there is 
not necessarily a galaxy within the NFW scale radius R s 
of the X-ray centroid, so MMGG sca i e and BGG sca i c are 
not necessarily defined for all galaxy groups. However, 
in the case of our clean sample, each group has at least 
one member within this radius so we do not have to deal 
with undefined centers. Uncertainties in the galaxy po- 
sitions are much smaller than the sizes of the galaxies, 
and are therefore negligible compared to the offsets from 
halo centers that we are capable of measuring with weak 
lensing. 

3.2. Centroid Candidate Centers 

The central galaxy is not always observationally ob- 
vious, and selection of an incorrect galaxy can produce 
statistically undesirable results when studying a sample 
of groups or clusters. This problem has motivated the 
use of centroids based on the positions of some or all 
group members, which can be weighted by their prop- 
erties such as flux or stellar mass, with the hope that a 
robust statistic can be less pr one to large errors than the 
choice of a single galaxy (e.g.,|White et al.|1999[ Carlberg 



et al. 2001 Berlind et al. 20061 IJee et al.||2011). 

J I ; II 1 1 II r 
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Additionally, ot her probes of groups and clu sters such 
as X-ray and SZ (Sunyaev & Zeldovich 1972) observa- 
tions of hot gas, and gravitational lensing can be used to 
find halo centers. Deep pointed observations can map the 
gas distribution in great detail for bright systems that 
are nearby or massive, but centering uncert aint ies can 
be significant for fainter systems (see Section 2.1). Simi- 
larly, only very massive systems produce a large enough 
lensing signal to study their spatial mass distribution in- 
dividually, and lower mass systems (like the groups stud- 
ied here) require stacking, such that centroids cannot be 
determined for each individual group from lensing alone. 

Here we consider four centroid candidates: 

• CN: the centroid of member galaxies. 

• CM: the centroid of member galaxies weighted by 
stellar mass. 

• CF: the centroid of member galaxies weighted by 
flux 

• X-ray: the X-ray centroid. 

Uncer tain ties on the X-ray positions were discussed in 
Section 2.1 and have a mean value of 22" or 134 kpc. 
For the other centroid candidates (CN, CM, CF), the 
coordinates are computed using a weighted mean, 



i=l 

N ' 



(4) 



where Xi is the pair of coordinates (R.A., Dec.)i for each 
galaxy i, N is the number of group members, and Wi 
is the appropriate galaxy weight for each center defi- 
nition; Wi = l,M* j, /j for CN, CM, CF, respectively, 
where M*^ is the stellar mass and fa oc i[p - 4mi with 
fi and rrii the flux and apparent F814W magnitude for 
each galaxy. We estimate the errors for these weighted 
means using bootstrap resampling from the list of mem- 
ber galaxies. With an average of 26 member galaxies 
per group, the mean statistical uncertainties on the pro- 
jected galaxy centroid positions are 45, 52, and 50 kpc 
for candidates CN, CM, and CF, respectively, where we 
have taken the geometric mean of the uncertainties in 
two dimensions and converted the angular uncertainty 
to a projected physical distance at the redshift of each 
group. Groups with a higher projected density of mem- 
ber galaxies tend to have smaller centroid uncertainties, 
but improvements appear to be limited by contamina- 
tion in the outskirts from correlated structure (see Pa- 
per I for tests of purity and completeness). We have 
tested different centroiding schemes including iterating 
until the centroid and member list converge or restricting 
to red galaxies, but achieved qualitatively similar results 
as with the centroids presented. 

3.3. Offsets Between Candidate Centers 

Our first test of these various centers is to see how 
well they agree with one another. Figure [T] shows the 
distribution over all groups of the angular and physi- 
cal distance offsets between pairs of candidate centers, 



along with the distribution of uncertainties in centroid 
positions. Immediately we see that candidate centers do 
not always agree. For example, in 22% of groups the 
brightest galaxy within i?200c is not the most massive 
galaxy and the candidates are separated by up to several 
hundred kpc. The agreement level among pairs of galaxy 
candidates is typically 70 — 80% with a long tail in the 
distribution extending out roughly to the virial radius for 
these groups. The galaxy candidates are typically offset 
from the centroid candidates by 50— 100 kpc, again with 
tails of a few hundred kpc, and the centroid candidates 
are in slightly better agreement among themselves. The 
offsets between the X-ray centroid and other candidate 
centers are generally consistent with the statistical uncer- 
tainties on the X-ray centroid. When comparing galaxy 
centroids (CN, CM, and CF) to other candidate centers, 
the typical offsets are roughly consistent with the mean 
uncertainty on the centroid position, but there are long 
tails in the offset distribution that exceed typical uncer- 
tainties. 

These results are generally consistent with offsets 
found in other groups and cluster samples, though di- 
rect comparison is difficult given the variety of methods 
and data used for identifying objects and their centers. 
For ~ 30% of optical ly-selected groups in a similar mass 
range as our sample, Skibba et al. (2011 1 found that the 
brightest galaxy was not the central one, based on the 
relative positions and velocities of other member galax- 
ies. This is comparable to the level of disagreement we 
find between our galaxy candidates, for which choosing a 
central galaxy is ambiguous. Comparing the positions of 
BCGs to X-ray centr oids in 42 optically-selected clusters, 
Sheldon et al. ( 2001 ) noted a mean offset of 85 kpc. With 
an expande d sample of 94 clusters with matching X-ray 
detections, Koester et al. (2007a) found a very similar 
median BCG-X-ray offset ol 81 kpc, and noted ~ 20% of 
systems with offsets of several hundred kpc which were 
mostly due to confusion in identifying the X-ray position 
or BCG. In a stud y of 65 massive cluster s with higher 
quality X-ray data, Sanderson et al.| ( |2009[ ) found BCG- 
X-ray offsets of typically less than a few tens of kpc, with 
a few outliers that were merging systems. Our X-ray 
offsets are somewhat larger due to the statistical uncer- 
tainties in the centroid positions, with a mean (median) 
offset of 104 kpc (78 kpc) between the MMGG sca i c and 
X-ray centroid. 

4. WEAK LENSING METHODOLOGY 

4.1. The Approach 

Our stacked weak lensing approach to test candidate 
centers is sketched in Figure [2] The left column shows 
two separate galaxy groups (red ellipses) and their corre- 
sponding shear maps measured from the shapes of back- 
ground galaxies (gray sticks). For each group, two can- 
didate centers are defined (blue triangle and purple di- 
amond), and the shear maps are stacked around these 
positions in the middle column. The rightmost panels 
show the resulting lensing signal AS as measured radi- 
ally from the candidate center. We emphasize that the 
lensing signal for individual groups studied i n this paper 
is noisy (si gnal-to-noise ~ 1; see Figure 1 of Leauthaud 
et al.|[2010 ) , so we cannot directly identify the centers of 
weak lensing maps for individual groups and must stack 
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Fig. 1. — Distribution of projected offsets between pairs of candidate centers in our group sample, measured in arcseconds (upper right; 
red) and Mpc(lower left; blue). The angular and physical offset distributions are not identical because the groups span a range of redshifts. 
The fill ed p urple histograms on the diagonal p anels show the distribution of statistical uncertainties for each centroid position, described in 
Section |2.1| for the X-ray centroid and Section |3.1| for the others. The y-axis gives the fraction of groups in each bin; bin sizes are 50 Mpc 
(bottom left and diagonal) and 10" (upper right). 

many groups; in this sense Figure [2] is exaggerated. 

Qualitatively, the amplitude oithe lensing signal is 
maximized when the lens position used for stacking co- 
incides with the true center of mass in each system, and 
the signal is suppressed when the nominal position devi- 
ates from the true center of mass. The two curves agree 
at radii much larger than the typical centering offset. 
More formally, the lensing signal around miscentered ha- 
los was first studied in the context of satellite galaxies 
dNatara.jan fc Kneib| [l997l JH ~' 1 ! ' 

|Seijak|irf)Oy||Yang et al.|iJ(J0iJ 



H udson et aL||1998| [Guzik fc 
2U06]), and later applied to 



the prob lem of uncertain group centers ((Johnston et al. 
2007a|b[ ). 



(J bservationally, our aim is to find the candidate center 
that maximizes the lensing signal on small scales, indicat- 
ing that it best traces the center of mass. Furthermore, 
we would like to model the signal to infer the underlying 



mass profile and the typical offsets between our tracers 
and the true center of mass in halos. Interpreting the sig- 
nal is complicated by a number of effects including the 
shapes of halo profiles and the properties of galaxies and 
subhalos which will be discussed further in Section [6] 

4.2. Models 

To interpret the mean lensing signal AS defined in 
Equations (T]) and p]), we construct a model for the sur- 
face mass density zJ[R) of a typical lens. Contributions 
to in our group sample come primarily from the 

dark matter halo and the central galaxy (if the center is 
defined to be at the position of a galaxy). The parame- 
ters of the model are introduced in Table [T] We include 
three variants of the model: a centered version where 
the halo center is fixed at the position of the candidate 
center, an offset model which allows a distribution of off- 
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Fig. 2. — Schematic illustration of stacked lensing around different candidate centers. Candidate centers are denned in each group (left), 
then shear maps are stacked around each position (middle), and azimuthally averaged to compute AS profiles (right). 



sets between the candidate and true halo center, and the 
full model which adds freedom to the halo profile and 
allows for excess mass in the form of a subhalo around 
the candidate center. 

We model the average mass density in halos with a 
spherical NFW profile, for which the project ed surface 
density Ehaio(-R ) = ^nfw(R) is given in e.g. Wright & 



Brainerd (2000), with halo mass and concentration as 
two free parameters. For the centered and offset models 
we will assume a mass-concentration relation from Zhao 



et al. ( |2009[ ), leaving mass as a single free parameter for 
the halo component, while both mass and concentration 
are free parameters in the full model. 

When the surface density of a spherically symmetric 
halo is measured around the correct center of mass, we 
have S(-R) = Shaio(-R)- If there is an offset i? D ff in the 
lens plane between the true center and the position used 
for measurement, the surface density measured at the 



coordinates (R, 9) relative to the offset position is ( Yang 
et al.||2003[ Appendix B) 



£fi o (i?,0|i? off ) = S halo i? 2 + R 2 oS - 2RR off cos^ . 

(5) 

The azimuthally-averaged surface density around the off- 
set position is 



1 f 27T 

S halo(^l^offO = / d^rflo^^off) 

27T Jo 



(6) 



For an e nsemble of halos with a di stribution of offsets 
P(Ros), Johnston et al. ( |2007a|b| generalized Equa- 
tion ([6]) to give the mean azimuthally-averaged surface 
mass profile stacked around the offset positions 



E halo (i?|P(i? off )) = 



-off 



P(Ros) Ehalo(-^l-ftoff) d-R, 



off- 



(7) 



The mean surface density inside a radius R is 



0< R\P(R oB )) ' 



R 2tt oo 



ttR 2 







P(Ros) ££ alo (i?', 6»|_R off ) 



x R'dR'dOdRos 



= ^2 1 z°L(R'\P(Ron))R'dR'. 

(8) 

To model the lensing signal from a large sample of 
galaxy clusters centered around BCGs, jJohnston et al.| 
(2007b) used a distribution of offsets P(-K ff) estimated 
from mock catalogs. In their model a fraction of BCGs 
correctly identified the centers of halos (i? ff = 0) , while 
the remaining clusters had a distribution of offsets given 
by 

P(R oS ) = ^exp(-ff). (9) 

(Toff V 2 °nff / 



This model, called a two-dimensional Gaussian or a 
Rayleigh distribution, was chosen based on mock cata- 
logs to which their cluster-finding algorithm had been ap- 
plied. The mocks suggested that the fraction of correctly- 
centered clusters depended on richness, with higher rich- 
ness clusters more likely to be centered correctly. In other 
clusters, the central galaxy was not correctly identified 
as the BCG, and the distribution of offsets between the 
BCG and true halo center could be described by Equa- 
tion ([9]) with the parameter a Q g — 420 h" 1 kpc describ- 
ing trie typical offset scale, independent of cluster rich- 
ness. 

We can think of the offset more generally in three di- 
mensions, where we assume the offset in each dimension 
is normally distributed with mean zero. The observed 
offset in the line-of-sight dimension might not have the 
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TABLE 1 
Model Parameters for AE(_R) 



Parameter 


Description 


Centered 


Offset 


Full 


Prior mean 


Prior a 


Restrictions 


log(M2OOc/M ) 


halo mass 


free 


free 


free 


13.5 


0.8 




C200c 


halo concentration 


tied 


tied 


free 


4.0 


.3.0 


1 < C200c < 10 


foff (kpc) 


offset distance scale 


fixed 


free 


free 





200 


CT ff > 


log(M gal /M Q ) 


stellar mass 


fixed 


fixed 


fixed 


<log(M*, gal /M )> 






log(M sub /M ) 


subhalo mass 


omitted 


omitted 


free 


(log(Af*, g ai/M )> 


1.0 


log(M sub /M ) > 10 



Note. — For each model, ga i is fixed to the mean photometrically-estimated stellar mass of the central galaxy (for 
galaxy candidates) and zero otherwise. For the centered and offset models, halo concentration is set by the relation of |Zhao| 
|et al.| p009) . 



same variance as the dimensions in the lens plane because 
of redshift-space distortions, but as long as the distribu- 
tion in three dimensions is joint-normal and the vari- 
ance in the two dimensions of the lens plane is equal, the 
projected offset distribution will take the form of Equa- 
tion @ after marginalizing over the line of sight. 

The miscen tering component employed by |Johnston 



et al. ( 2007b I added two free parameters to their model: 



the traction of miscentered groups, and the scale length 
of the offset distribution. These additional parameters 
could not be well-constrained by the data, and had to be 
constrained using strong priors from mock catalogs. 

Our approach is to use an offset model with a single free 
parameter to constrain the scale of the miscentering dis- 
tribution empirically. We do not include separate compo- 
nents for cent ered and miscentered groups as |Johnston"| 
et al. ( |2007b[ ) did, as the single parameter model gen- 
erally provides sufficient fits to the current data. Also, 
correctly chosen central galaxies may still be offset from 
their halo centers, an effect that was not included in the 



mocks of Johnston et al. ( 2007b I but has bee n tentatively 
seen in lensing maps of individual clusters ( |Oguri et al. 
20101. 



The second component of the model is the surface den- 
sity contributed by the central galaxy, £ ga i(i?). The 
shape of this mass profile is uncertain, but its contri- 
bution is subdominant even at the smallest radii that we 
probe (R ~ 50 kpc). For the centered and offset mod- 
els, we simply model the galaxy component as a point 
source, E ga i(i?) = M ga i/(nR 2 ), with M ga i fixed to the 
average stellar mass M*, ga i of the central galaxies in the 
ensemble as estimated from their SEDs. The centroid 
candidates (CN, CM, CF, and X-ray) typically do not 
have a galaxy very near the center, so we do not include 
a contribution from E ga i when modeling the signal from 
these candidates. 

We can write the centered and offset models as the sum 
of the halo and galaxy components: 



AS cen 
AE' 



(R): 
° S (R): 



:A£ NFW (i?) + A£ ga i(i?) (10) 
:AS^ w (i?|P(i? off )) + AS gal (i?). (11) 



Note that AS ccn (i?) = AE off (R\R oB = 0). We later test 
a scenario with the full model in which the candidate cen- 
ter has an additional mass component in the form of a 
dark matter subhalo that is more extended than the stel- 
lar profile which will be described further in Section [5.3| 
Though the measured signal is an ensemble average 
coming from halos with a range of masses and mass pro- 
files, our model consists of a single profile for simplicity 
and because the range of halo masses inferred from group 



X-ray luminosities is relatively small. We restrict our 
analysis of the lensing signal to R < 1 Mpc, where the 
effects of halo truncation and correlated structure should 
be unimportant to the lensing signal, and note that this 
adequately covers the range of centering offsets shown 
in Figure [TJ At small scale s, the assumption of weak 



shear becomes less accurate. Mandelbaum et aLI ( |2006a I 
derive a correction term to the surface density contrast 
that depends on E crit and E(i?) for the sample. These 
correction factors have b een computed for this sample by 
Leauthaud et al. (2010) and shown to be fairly small on 
the scales we probe (of order 10% of the measured signal 
in our innermost bin for a good center), so we do not 
include them in our analysis. 

To fit the models to the data, we attempt to find the 
parameters that minimize 



= E 



(AEdata(i?i) 



my 



(12) 



where Oi is the measurement uncertainty on AEdata(-Rj)- 
In practice, we use a Markov Chain Monte Carlo 
(MCMC) approach to efficiently explore the multi- 
dimensional parameter space maximizing the logarithm 
of the likelihood, L oc exp(— x 2 /2). We employ Gaus- 
sian or lognormal priors for each parameter, with means 
and standard deviations given in Table [T] The data are 
unable to constrain a lower limit to the subhalo mass or 
an upper limit to the concentration, so we restrict these 
parameters to log(M su b/M ) > 10 and C200C < 10 when 
they are free. We also require C200C > 1 and er ff > 0. 

5. RESULTS 
5.1. Weak Lensing on the Full Sample 

Given that Figure [T] shows a wide range of offsets be- 
tween different choices of group centers, we proceed to 
test how well the candidate centers trace the underly- 
ing matter distribution. We begin by studying the full 
sample of groups in this section, and in the next section 
we focus on the subset of groups with significant offsets 
between candidates. 

For each center candidate, we measure the lensing sig- 
nal AE(i?) from all groups in annular bins around the 
center. The results are shown for the eight candidate cen- 
ters in Figure JS] Each panel represents a different can- 
didate center from Section[3j for which we plot the mea- 
sured AE(iZ) (black pointsjalong with models AE ccn (i?) 
(thick blue) and AE off (i?) (thin magenta). For illus- 
tration, we show the halo and galaxy components of 
the models for MMGG sca i : AEnfw(-R) (green dashed), 
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AE^w(-Rl-Roff) (orange dot-dashed), and AE gal (i?) (red 
dotted) . 

The mean and standard deviation (la) of the pos- 
terior distribution for each parameter, along with the 
mean central galaxy masses and \ 2 values for the fits, 
are shown in Tabled The x 2 values have not been nor- 
malized by the number of degrees of freedom v. Each 
model is fit to 6 data points with 1 and 2 free parame- 
ters for the centered and offset models, respectively, so 
v = 5 for the centered model and v — 4 for the offset 
model. When \ 2 w v, the model is consistent with the 
data given the uncertainties. A high value of x 2 indi- 
cates that the model does not fit the data; for example, 
the data are inconsistent with the model at a 95% confi- 
dence level when x 2 > 11.07 (9.49) for v = 5 (4). 

There are clear differences among the lensing signals 
for the eight candidate centers. The signal generally rises 
toward small scales for the galaxy candidates, while there 
is a turnover for the centroid candidates indicating that 
they are poorer at tracing the center of mass in these 
groups. The stellar mass in the candidate central galax- 
ies produces some of the difference in the signal, but 
the best-fit models show that the galaxy candidates have 
smaller offsets from the halo centers when compared with 
the centroid candidates even when accounting for the 
stellar mass. We will explore the possibility of ad ditio nal 
mass in the central galaxy candidates in Section 5.3 It 



is important to note that the offsets measured from the 
lensing signals for the X-ray position and other centroid 
candidates are primarily due to the large statistical un- 
certainties in the centroid positions, shown in Figure [T] 
Figure [3] shows the impact of these centroid uncertain- 
tics on the lensing signal, but wc cannot infer from this 
data that there are significant intrinsic offsets between 
the true centroid positions and the center of the dark 
matter halo. 

For each of the galaxy candidates, the x 2 value in Ta- 
ble [2] indicates that the centered model provides a good 
description of the data given the error bars. The extra 
degree of freedom in the offset model does not signifi- 
cantly improve the fit, and the constraints on a Q s are 
consistent with zero at the 2a level. The best-fit values 
for a Q ff and the quality of the fits are statistically con- 
sistent for each of the galaxy candidates, so we do not 
strongly favor one candidate over another, though there 
is marginal evidence that the candidates defined within 
an NFW scale radius of the X-ray position (MMGG sca i 
and BGG sca io) trace the halo center slightly better than 
candidates defined over the larger area within i?200c- 

On the other hand, the centered model generally does 
not provide good fits to the data for the centroid candi- 
dates, and the offset models do improve the fits. The 
best-fit offset parameter a s is significantly larger for 
the centroid candidates than for the galaxy candidates 
and deviates from zero by more than 2a for three of the 
centroid candidates. The offset model fits two centroid 
candidates (CM and CF) only marginally well so these 
would perhaps be better fit by a more complicated offset 
distribution. 

The effect of miscentering on the halo mass constraints 
is comparable to the statistical uncertainty for this sam- 
ple. The best-fit halo masses assuming a centered model 
for the centroid candidates tend to be ~ la lower than for 
galaxy candidates. Halo masses increase slightly when 



offsets are allowed in the model, and candidates with 
larger offsets are more affected. When miscentering is 
not accounted for in the model, halo masses are under- 
estimated by 5 — 30% compared to when we include the 
effect, depending on the choice of center. 

We have also tested the effect of fixing the concen- 
tration parameter with our assumed mass-concentration 
relation, and freeing this parameter does not qualita- 
tively change the results. We do not obtain a good con- 
straint on the concentration, but the best-fit value for the 
galaxy candidates is consistent with our assumed value. 
The centroid candidates prefer somewhat lower concen- 
trations, but the values are still consistent within the 
large error bars. 

5.2. Groups with Discrepant Candidate Centers 

While the stacked lensing signal for the full sample 
of groups shown in Figur e [3] exposes differences between 
candidate centers, Figure [1J shows that the candidate cen- 
ters are identical in many groups. In order to more di- 
rectly compare different centering schemes, wc now select 
only the groups in which candidate centers disagree by a 
measurable amount. Figures [4] and [5] present this analy- 
sis, with each row showing the lensing signal centered on 
the MMGG S caie (left column) and a different candidate 
(middle column) for the subset of groups where the two 
candidate centers are separated by more than 50 kpc. 
The rightmost column shows the histogram of offsets be- 
tween the two centers repeating Figure [l] on logarithmic 
axes. The number of groups and the mean redshift of 
each subset are shown at the top right. We have chosen 
the MMGG sca i as a fiducial center here since Table [2] 
suggests that it is among the best candidates, but using 
the BGG sca ic gives a consistent picture. 

In each row of the figure, we first model the lensing sig- 
nal in the lefthand panel assuming that the MMGG sca i c 
correctly traces the halo center, fitting for the halo mass 
with concentration fixed. Next, we take that best-fit cen- 
tered model and convolve it with the distribution of off- 
sets in the righthand panel to model the signal in the 
middle panel. This is equivalent to replacing the offset 
distribution of Equation ([9]) with the empirical distribu- 
tion of offsets between two candidates, leaving no free 
parameters for the models shown in the middle panels. 
We again add a central point mass fixed to the mean pho- 
tometric stellar mass of the candidate galaxies for each 
of the galaxy centers. Though we expect that galaxies 
offset from the halo center have subhalos, the lensing 
data are currently insufficient to constrain such an addi- 
tional contribution. The halo mass is fitted using only 
the lensing data centered on the MMGG sca i c , so the level 
of agreement between the offset model and the data in 
the middle panel can be read as a consistency test. 

Looking first at the comparisons between the 
MMGG SC aie and the centroid candidates in Figure |4j we 
see that the lensing signal is reasonably well-described 
with an NFW halo centered on the MMGG sca ie, with a 
mean halo mass consistent with that measured for the 
full sample in Table [2] Additionally, the offset model 
gives a decent representation of the lensing signal mea- 
sured around the centroid candidates shown in the mid- 
dle panels. We interpret these results as indicating that 
centroid candidates do not trace halo centers as well as 
galaxy candidates. 



10 



George et al. 



o 

Q- 

© 



DC 

W 

rr 

IW 



< 




0.1 1.0 0.1 

Physical transverse distance, R (h 72 



1.0 

Mpc) 



Fig. 3. — Weak lensing signal stacked on the full sample of groups around different centers, along with centered (thick blue, AE con ) and 
offset (thin magenta, AE off ) models. Halo and central components of these models are shown for MMGG sca i e (green dashed for ASJj™ w , 
orange dot-dashed for AS^, W , red dotted for AS ga i). The top row shows the signal around galaxy candidates, while the bottom row shows 



ccntroid candidates. The signal (black points) is measured in radial bins with the first spanning 20 
then logarithmically spaced from 70 kpc to 1 Mpc. 



70 kpc for sufficient signal-to-noise, 



TABLE 2 

Parameters Constraints for AS (J?) 



Centered Model 



Offset Model 



Candidate (log(M* iga i/M Q )) 


log(M 2O oc/M ) 


C200c 


X 2 (v = 5) 


log(M 2OO c/M ) 


C200c 


<r off (kpc) 


X 2 {v = 4) 


Galaxy candidates 
















MMGG sca i c H.3 


13.43 ± 0.07 


4.0 


4.2 


13.45 ± 0.07 


3.9 


18.2 ± 11.3 


4.2 


MMGG R20 o 11-4 


13.43 ± 0.06 


4.0 


5.3 


13.47 ± 0.07 


3.9 


28.3 ± 15.2 


4.4 


BGG sca le 11.2 


13.40 ± 0.06 


4.0 


2.9 


13.43 ± 0.07 


1.0 


16.8 ± 10.0 


2.8 


BGGr,200 11-3 


13.42 ± 0.07 


4.0 


4.4 


13.45 ± 0.06 


3.9 


24.8 ± 12.0 


3.0 


Centroid candidates 
















CN 


13.35 ± 0.07 


4.0 


21.0 


13.51 ± 0.07 


3.9 


65.0 ± 18.0 


2.8 


CM 


13.39 ± 0.07 


4.0 


10.5 


13.45 ± 0.08 


3.9 


34.5 ± 21.6 


9.3 


CF 


13.38 ± 0.07 


4.0 


14.0 


13.51 ± 0.09 


3.9 


66.9 ± 29.9 


9.5 


X-ray 


13.30 ± 0.08 


4.1 


15.2 


13.42 ± 0.08 


4.0 


57.1 ± 16.9 


3.5 



Note. — For both models, M* ga j is fixed to the stellar mass of the central galaxy (for galaxy candidates) and zero otherwise, 
and the halo concentration is fixed by the relation of |Zhao et al.| l [2009| . 



In the samples where the MMGG sca i c disagrees with 
other galaxy candidates (Figure [5j), the lensing signal is 
different than in the centroid comparison. The best-fit 
halo masses are significantly lower than in the full sample 
measured earlier. The lensing signals are noisier, in part 
because of the smaller sample sizes, and one case deviates 
significantly from the offset model. Though it appears 
that the BGG sca io produces a higher lensing signal than 
the MMGG sca io at small radii in this direct comparison, 
that profile is not well-fit by a centered NFW model and 
the fitted mass is similarly low. 

Figure [6] shows how two samples of groups with dif- 
fering galaxy candidates are distributed in redshift and 
X-ray luminosity relative to the full sample of groups. 
Cases where the most massive group member lies in 



the outskirts (MMGG sca i ^ MMGG R2 oo) appear to be 
evenly distributed throughout the sample. Groups where 
the brightest galaxy near the X-ray position is not the 
most massive (MMGG sca i c 7^ BGG sca i ) tend to be at 
higher redshifts, which can also be seen in the mean red- 
shifts for the samples in Figure [5j This illustrates how 
an observed-frame selection of BCGs, like that used for 
BGG sca ie here, tends to pick up bluer star-forming galax- 
ies at higher redshifts. 

The low masses shown in Figure [5] could be attributed 
to statistical fluctuations with the small sample size 
(24 — 40 groups, as compared with 91 — 98 groups in 
Figure H|. To test this idea, we performed jackknifc 
tests using randomly selected samples of the same num- 
ber of groups without replacement. In each of 1000 
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Fig. 4. — Lensing signal for subsets of groups where centroid candidates are offset from MMGG sca i by > 50 kpc. The left column shows 
the signal stacked around MMGG sca i e , the middle column shows the signal stacked around an alternate candidate, and the right column 
shows the projected distribution of offsets between the two candidate positions. Models are discussed in the text with line styles and colors 
as in Figure and fit parameters stated within the left column. The number of groups with large offsets used in each row is stated in the 
right column, and the gray hashed boxes indicate the number of groups excluded from this analysis because the two candidates agree to 
within 50 kpc. 



random samples we measured the lensing signal around 
MMGG sca i e and fit a model with a central point mass 
equal to the average stellar mass and a centered NFW 
with halo mass as a free parameter and concentration 
fixed. Of the 1000 random samples of 24 groups, 19% 
had a best-fit halo mass lower than that found for the 
sample of 24 groups where MMGG scak! ^ MMGG R20 o- 
In samples of 25 groups, none had a lower best-fit halo 
mass than log(M 2 oo c /Mq) = 12.57, the value found for 
the sample with MMGG sca io 7^ BGG sca i . And in sam- 
ples of 40 groups, only 1% had a best-fit mass lower than 
the sample where MMGG sca i 7^ BGGr200- For a given 
X-ray luminosity, groups with discrepant galaxy candi- 
date centers appear to have lower masses than the rest 
of the sample and in some cases more disturbed mass 
profiles. 

Figure [6] shows a redshift dependence in the fraction 
of groups where the brightest galaxy (measured from the 
observer- frame F814W flux) does not have the highest 
stellar mass. High redshift groups may have more dis- 
turbed lensing profiles for different reasons than hav- 
ing ambiguous galaxy center candidates, so we have re- 
peated the jackknife tests restricting to groups at z > 0.7. 
Even among high redshift systems, the groups with a 
brighter galaxy than MMGG sca i e are outliers; less than 
1% of random samples of groups produce a best-fit mass 



lower than the samples where MMGG sca i c 7^ BGG sca i 
or MMGG sca i c 7^ BGGR200, while just 17% of random 
samples are fit by lower masses than the groups with 

MMGGscale 7^ MMGGr,200- 

Visual inspections of the groups with discrepant galaxy 
candidates have not revealed obvious differences from the 
rest of the sample. Given the noisy lensing signal mea- 
sured around such groups, one might worry that they are 
not real associations but chance projections. However, 
these groups are not preferentially found near the X-ray 
flux limit shown in Figure [6] and they do not have un- 
usually low numbers of member galaxies associated with 
them, so they are unlikely to be impurities in the group 
sample. 

Groups with multiple massive galaxies and disturbed 
density profiles are likely to have merged recently, a point 
we discuss further in Section [U 

5.3. Model Uncertainties 

Having established that the lensing data can be used to 
constrain the choice of a tracer for the centers of mass of 
halos, we would like to know how well a given candidate 
actually traces halo centers. In this section we explore 
a more general model to see how certain assumptions, 
namely the mass-concentration relation and the form of 
the central mass component, affect our results. 
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Fig. 5. — Lensing signal for subsets of groups where galaxy candidates are offset from MMGG sca i e by > 50 kpc. Plot columns and style 
are the same as in Figure [4] 



44.0 



43.5 



w 43.0 

O) 

Si, 

"*42.5 



O 



42.0 

41.5 




• • 

♦ 


> 

• ♦ : • • * • - ■ 






. ; ■ ■ 






^ MMGG sca le ^ MMGGr200 


• ✓ 

* 

* 

* 


. MMGG scalc / BGG scaJe 



be a satellite or the system may be unrelaxed. Detailed 
modeling of strong and weak lensing observations of in- 
dividual clusters has found evidence for subh alos around 
member galaxies even near cluster centers ( jNatarajan 
et al. 2007 [2009). Additional mass in the central term 



^gai(^j could hide significant offsets between the galaxy 
candidates and host halos, effectively filling in the decre- 
ment at small radii in the halo term Ehaio(-R)- 

To address this degeneracy, we apply our full model 
with an additional mass component around the galaxy 
candidate, in addition to the point source representing 
the stellar mass. We use the functional form of a trun- 
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Fig. 6. — Distribution of redshifts and X-ray luminosities for the 
group sample. Cases where different galaxy candidates disagree 
are identified according to the legend. The dashed curve shows the 
4a X-ray flux limit of 1.0 X 10 5 ergcm - 2 s — 1 reached in 96% of 
the field converted to a limiting luminosity. 

The results shown in the previous section suggest that 
galaxy candidates trace the halo centers better than cen- 
troid candidates, because the stacked lensing signal on 
small scales is greater when centered on a galaxy. While 
the model used to fit the signal for galaxy candidates in- 
cludes a component to account for the stellar mass that 
exists in galaxy candidates but not in centroid candi- 
dates, it is plausible that we have underestimated the 
amount of mass in the central galaxies, either in stars 
and baryons or in a dark matter subhalo. In a typi- 
cal halo model the central galaxy has no subhalo, so in 
the latter scenario the galaxy candidate could actually 



cated non-singular isothermal sphere (Pastor Mira et al. 
2011] ) 

P^"> = I -| I 2/ 2 i 2/2 \ ' ( 13 ) 

We fix the core radius r corc = 0.1 kpc and set the trun- 
cation radius r cut to the distance at which the local mass 
densities of the subhalo and halo components are equal, 
as measured along the line connecting their centers as- 
suming an offset distance er ff . To facilitate comparisons 
with the stellar mass, we cast the normalization of the 
model in terms of a free parameter M su b (instead of po) 
which we will refer to as the subhalo mass, defined as the 
mass enclosed within 5 kpc, the mean half-light radius for 
our sample. Since the core radius is much smaller than 
our innermost lensing measurement, the exact value cho- 
sen does not influence our results. The choice of trun- 
cation radius and functional form of the subhalo mass 
profile does affect the mass normalization, but the qual- 
itative shapes of the parameter degeneracies are consis- 
tent for a variety of truncation schemes and for mass 
profiles that are isothermal, or nearly so, across the in- 
ner few tens of kpc. However if the subhalo component 
is not truncated at all, or if it is mo deled instead as an 
NFW profile as Pastor Mira et al. find in simulations 
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with subhalos, its mass can become degenerate with the 
full halo mass. 

For a given lensing signal, a larger halo concentration 
can also compensate for some miscentering effects, so we 
free the C200C parameter to study the degeneracy with 
<7 ff without the restrictions of the model for concentra- 
tion used before. In practice, concentration is not well- 
constrained by this data set, so a prior is still needed to 
restrict the range of values when fitting. 

An additional model uncertainty comes from the form 
of the distribution of offsets P(R a tf) that is used. The 
offsets between galaxy cand idates in Figure [T] appe ar to 
support the model s used by Johnston et al. (2007b) and 
Oguri et al. (2010), with a significant fraction of agree- 
ment between candidates along with a wide tail in the 
distribution. To avoid the need for additional free pa- 
rameters describing both well-centered and miscentered 
populations in the offset distribution, we attempt to se- 
lect a clean sample of groups and constrain the scale of 
offsets between the candidate galaxy center and the halo 
center assuming the single distribution given in Equa- 
tion (p)]). We test this model on the sample of groups 
for which the four galaxy candidates agree, i.e. where 
the brightest and most massive group members are one 
and the same, and where this galaxy lies within an NFW 
scale radius of the X-ray centroid. There are 85 groups 
with such unambiguous galaxy centers. 

The lensing signal for these groups stacked around the 
galaxy candidate is shown in Figure [7] To summarize our 
model, we fit the signal with an NFW halo and a trun- 
cated isothermal subhalo of the form in Equation ( 13 ) in 



addition to a point mass fixed by the photometric stel- 
lar mass. The point mass is fixed at the center of the 
subhalo, which is allowed to be offset from the center 
of the halo. This leaves four free parameters: the halo 
mass Af 2 ooc and concentration C2000 the offset scale 0~ o ff 
between the galaxy and halo, and the subhalo mass nor- 
malization Af su b- 



250 




0.1 1 

Physical transverse distance, R (h^ Mpc) 

Fig. 7. — Stacked lensing signal for 85 groups where galaxy can- 
didates agree. Model curves show the maximum likelihood solu- 
tion for an offset NFW halo (orange dot-dashed), subhalo (purple 
dashed), and a point source for the stellar mass (red dotted), along 
with the sum (thick solid blue). There are degeneracies among 
some fit parameters; for example, a model with smaller centering 
offsets and a less massive subhalo also fits the data. 

Figure [8] shows the results of the MCMC analysis ex- 



ploring this parameter space. Each panel with blue con- 
tours shows the 68% and 95% regions of the joint pos- 
terior probabilities for a pair of parameters, marginaliz- 
ing over the other parameters. The top panel in each 
column shows for each parameter the arbitrarily nor- 
malized prior from Table [l] (dashed green curves) and 
the one-dimensional posterior probabilities (black his- 
tograms) while marginalizing over the other parameters. 
Note that the data are unable to constrain a lower limit 
on the subhalo mass or an upper limit to the concentra- 
tion, so the posterior distributions are sensitive to our 
priors in those regions. 

The data are consistent with a model in which there 
are no subhalos around the candidate central galaxies, 
though a value comparable to the average stellar mass 
of the galaxies is preferred, and significantly more mas- 
sive subhalos are ruled out. We emphasize the stellar 
mass is already included in the model as a point source, 
so the subhalo component represents mass in excess of 
the stellar mass that was estimated from the photome- 
try. Recall that the M su b is defined as the subhalo mass 
within the t ypic al half-light radius of 5 kpc; integration 
of Equation ( 13 1 out to large radii for the maximum like- 



lihood parameters gives a total subhalo mass of roughly 
five times the observed stellar mass component. The 
truncation radius for this model is 21 kpc. The pref- 
erence for a nonzero subhalo mass may simply indicate 
that the photometric stellar mass is underestimateoJ^J 
or that the point source model for the stellar term is in- 
sufficient. A subset of correctly-centered groups (or an 
excess at small scales in the offset distribution relative to 
the model) could also explain the signal which we have 
fitted as a subhalo component. The halo mass is well- 
constrained and consistent with the values measured for 
the full group sample with the simpler model in Sec- 
tion [5]2_Thedlata prefer a higher value of C200C than the 
Zhao et al. ( 2009 ) prediction for this halo mass and red- 
shift, but the concentration is not well-constrained. For 
the offset distribution, the data are consistent with no off- 
sets between the central galaxy candidates and the halo 
centers, though a value of roughly 50 kpc is preferred. 
Large offsets (a > 100 kpc) are ruled out for an offset 
model of a single two-dimensional gaussian distribution, 
though individual outliers may exist. On average, mas- 
sive galaxies trace halo centers quite well. 

There are degeneracies between the offset distance 
scale and the subhalo mass as well as the halo concen- 
tration. If the subhalo mass within 5 kpc is smaller than 
~ 10 11 Mq, the offset scale must be smaller than 50 kpc, 
but for larger subhalo masses the offset scale can reach 
100 kpc. A constraint on the total mass enclosed within 
small scales, perhaps from strong lensing or measure- 
ments of stellar velocity dispersions for the central galax- 
ies, could improve constraints on a g. Similarly, C200C 
and Coff are correlated, with higher concentrations cor- 
responding to larger offsets. Thus when modeling a given 
lensing signal, if miscentering is not taken into account 
it can lead to an underestimate of the concentration. 

6. DISCUSSION 

12 Stellar mass es used in this work have been derived assuming 
alChabrierl (12003} initial mass function (IMF). Recent suggestions 
or a steeper subso lar IMF for early type galaxies could explain thi s 
discrepancy (e.g., Auger et al. 2 010||van" Dokkum fc Conroy|20"l0) . 
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Fig. 8. — Posteriors for the four parameters of a general offset model discussed in the text, applied to the lensing signal from Figure J7] 
Blue contours show 68% and 95% regions for the pair of parameters noted along the axes, marginalized over the other parameters. White 
crosses mark the maximum likelihood parameters. Top panels show the posterior distributions for single parameters while marginalizing 
over the others (black histograms), along with arbitrarily normalized priors (green dashed curves). 



Finding the centers of mass for dark matter halos 
presents both an observational challenge and an oppor- 
tunity to study the interplay between halos and central 
galaxies. In this paper, we have presented a method 
to test different tracers of group and cluster centers by 
comparing the weak lensing signal stacked around their 
positions. Our approach can in principle be applied to 
any analysis that would be affected by miscentering. For 
instance, centering algorithms can also be tested using 
satellite kinematics by identifying the candidate that is 
nearest to the dynamical center when averaging over the 
ensemble. The spatial clustering of galaxies can also pro- 
vide centering information; in practice the peaks in the 
galaxy density field are generally used in defining optical 
cluster catalogs, but clustering data could also be used 
to optimize the determination of their centers. Addition- 
ally, these approaches to optimizing the centers of a halo 
catalog can be fed back into algorithms used to find halos 
and their centers, allowing a more probabilistic approach 
to deal with cases where centering is uncertain. 

With X-ray detected galaxy groups in the COSMOS 
field, we find that individual bright and massive galaxies 



trace the centers of halos better than the nominal X-ray 
centroid or the mean position of group members, even 
when weighted by luminosity or stellar mass. Offsets be- 
tween the X-ray centroids and candidate central galax- 
ies are primarily due to the large uncertainties in the 
X-ray positions, which can reach roughly 200 kpc. Cen- 
troids based on the mean positions of galaxies also have 
significant uncertainties, though some offsets exceed the 
estimated errors. More stringent constraints on intrin- 
sic offsets between centroid candidates and halo centers 
could be obtained with deeper, high-resolution X-ray or 
SZ data, or with a sample of massive clusters and more 
member galaxies. 

For each pair of galaxy candidates defined in Section [3j 
roughly 20 — 30% of groups have ambiguous centers. Ei- 
ther the most massive galaxy lies far from the X-ray 
position, or the brightest galaxy is not the most mas- 
sive. Only two-thirds (85 out of 129) of the groups show 
complete agreement between all four galaxy candidates, 
which underlines the importance of testing any choice of 
center in a group or cluster catalog. 

The results of Figure [5] and the jackknife tests de- 
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scribed in Section |5.2| suggest that groups with ambigu- 
ous galaxy centers have low masses for their X-ray lu- 
minosities or disturbed halo mass profiles. A connection 
between centering offsets and the dynamical sta te of ha- 



los has been seen in observational studies (e.g., Forman 
fe Jones||1982[ |Katayama et al. 2003; Sanderson et a 



2009[) and~simulations (e g.,|Cohn fc White||20051 |Poole| 
eFaL|[2006l |Macci6 et al.||20U7| |Meto et al.||2007| |Skibba| 
fc Macc i6|2011 ). The ability to identify unusual or unre- 



laxed groups and clusters with a simple observable such 
as having different candidate centers could prove use- 
ful when trying to select a clean sample of relaxed sys- 
tems. Further studies connecting lensing measurements 
with member galaxy properties such as the distribution 
of colors may provide a clearer picture of the impact of 
halo mergers on galaxies and star formation. Similarly, 
a more general analysis of the gaps in stellar mass or 
luminosity and the spatial separation between massive 
member galaxies would complement the restricted set of 
center candidates studied here. 

There are several reasons why candidate centers can be 
offset from halo centers. Interlopers or incompleteness 
in the group member sample can result in the wrong 
galaxies being selected. Satellite galaxies can be more 
massive or luminous than centrals due to scatter in the 
relation linking the observable property to the mass of 
a halo or subhalo. In Paper I, we presented tests of our 
group membership and centering algorithm on mock cat- 
alogs designed to reflect the real uncertainties in X-ray 
positions and photometric redshifts. In that analysis, 
77% of central galaxies were correctly identified as the 
MMGG sca i c . The most common failure mode was the 
selection of satellites as MMGG sca io because they had 
higher stellar masses than the true centrals, which hap- 
pened in 12% of cases. The frequency of this occurrence 
depends on the way we populate halos with mock galax- 
ies following the st ellar mass-halo mass relation of Leau- 
thaud et al. (2012), but we note that scatter in thi s rela- 
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remaining centering failures were evenly split between 
cases where X-ray errors misplaced the search region and 
photometric redshift errors scattered centrals out of their 
groups. 

Observational uncertainties, such as scatter in stellar 
mass estimates (~ 0.2 dex) and catastrophic errors in 
photometric redshifts (for < a few percent of objects in 
our COSMOS sample), add to this scatter and increase 
the chances that a satellite is incorrectly identified as 
the central galaxy. Furthermore, merging activity is not 
modeled in the mocks catalogs, so these effects together 
can account for the ~ 30% of observed groups where 
galaxy candidates disagree. The offsets seen between 
centroid candidates (CN, CM, CF, and X-ray) are consis- 
tent with being due primarily to the large observational 
uncertainties in their positions. 

The offsets measured in this paper are significantly 
smaller t han the distribution me asured from mock cat- 
alogs by Johnston et al. (2007b). This could b e du e in 
par t to our choice ot onset model; in Sections |5.1| and 
|5.3| we a ssumed a single Rayleigh distribution of offsets, 
whereas |Johnston et al.| separated groups that were cor- 
rectly centere d fro m those that were miscentered. Our 
aim in Section 5.1 was simply to estimate a typical offset 



scale for comparing candidate centers, and in Section [B~3] 
we selected a sample where centering seemed unambigu- 
ous to study up close the offsets of massive galaxies which 
could be s losh ing around halo centers. The offsets found 
in Section |5.3| are comparable to the smaller componen t 



of the offset distribution measured by |Oguri et al. ( 2010[) 



While they are fairly small, modeling degeneracies with 
the subhalo mass and halo concentration increase their 
uncertainty. Future analyses of halo properties may ben- 
efit from including a range of models for the distribution 
of offsets, if the data can constrain a larger number of 
parameters. 

Anot her important diff erence between our analysis and 
that of IJohnston et al.l is that the maxBCG clusters 
are detected as optical galaxy overdensities, whereas the 
COSMOS groups studied here are detected in X-ray 
emission, which traces the dense regions near the cen- 
ters of halos, offering a better starting point for finding 
the centers of halos. When comparing the effects of mis- 
centering for different cluster catalogs, it is also worth 
noting that more massive clusters are larger, so offsets of 
a given distance can be more easily detected than in less 
massive systems. 

Further investigations with simulations could improve 
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can be offsets between the position of the most bound 
particle, the mass density peak, and the centroid of mass 
with a given smoothing scale. Offsets between these po- 
sitions in simulations have been shown to correlate with 
the dynamical state of a halo, with larger offsets seen in 
less relaxed halos that have experienced a recent merger. 
Projecting the matter density in simulated halos to com- 
pare with the observed lensing signals could help explain 
the poor fits to NFW profiles and unusually low masses 
we obtain from samples with ambiguous centers. Simi- 
larly, a better understanding of the form and evolution 
of halo and subhalo mass profiles, due to effects like gas 
cooling or heating and tidal stripping, could improve our 
modeling. 

Finally, we consider the implications of misccntcring 
for cosmological analyses. The abundance of massive 
halos is sensitive to the amplitude of matter fluctuations 
and the growth history of the Universe, and the precise 
determination of grou p and cluster masses is a c ritical 
aspect of this probe. Mandelbaum et al. (20101 have 



studied the effects of miscentermg on mass estimates for 
massive halos using simulated data and analytical pro- 
files, assuming the distribution of o ffsets from the mock 
catalogs of |Johnston et al.| ( 2007b ) . For clusters at the 
upper end of our mass range (~ 10 M Q ), they find that 
the weak lensing mass is underestimated by 25 — 30% if 
miscentering is ignored, and that this effect is stronger 
for less massive halos (see their Figure 3). They also 
show that the accuracy of mass estimates depends on 
the assumed concentration as well as the inner and outer 
radius measured, and suggest excising the inner regions 
from the analy sis because of these uncer tainties. In a 
separate study, [Mandelbaum et al. ( 2008[ ) measured the 
mass-concentration relation tor several samples of galax- 
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ies and clusters, including the maxBCG sample. Miscen- 
tering of clusters is also a concern for measuring concen- 
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may be overestimated, based on the concentrations they 
derived from lensing and by comparison with the distri- 
bution of offsets betw een BCGs and X-ray peaks in a 



subsample of clusters (Koester et al. 2007a). Still, they 



had to trade off statistical power by excluding data at 
small radii to reduce systematic uncertainties from mod- 
eling the distribution of offsets. 

In this paper, we have endeavored to improve the ac- 
curacy of finding halo centers and to provide constraints 
on the distribution of offsets between observational trac- 
ers and underlying mass centers. With more accurate 
centers (and uncertainties on those positions), we can 
more reliably use data at small radii and improve statis- 
tical constraints from group and cluster surveys. Table [2] 
suggests that the halo mass inferred from the lensing 
signal stacked around our galaxy candidates is affected 
at the 5 — 10% level if miscentering is ign ored. This is 
significantly smaller than the bias seen by |Mandelbaum| 
et al. (2010) despite the trend that masses are increas- 
ingly underestimated for lower mass halos. We attribute 
the difference to the smaller offsets seen in this s ample 
compared to the m iscentering distribution used by |John- 



ston et al. (2007b) where a fraction of groups had a dis 
tribution of offsets with a s = 420ft -1 kpc. However, 
we also use a different model for the offset distribution 
and can see in Table [2] that the masses estimated from 
centroid candidates are more biased (15 — 30%) when 
miscentering is not addressed in the model. The statis- 
tical uncertainties on our halo mass estimates are still 
comparable to the centering bias, and with this sample 
we are currently unable to put a significant constraint on 
the halo concentration. But with larger group and clus- 
ter samples, our approaches to optimizing the choice of 
halo centers and modeling the distribution of offsets will 
enable better constraints on halo mass profiles for both 
astrophysical and cosmological applications. 

7. SUMMARY AND CONCLUSIONS 
We summarize the main results of this paper as follows: 

1. In our data set, different definitions of group cen- 
ters do not always agree and occasionally show 
large offsets from one another. Candidate centers 
based on the locations of the brightest or most mas- 
sive galaxy differ in 20-30% of cases with a wide 
range of offsets. Centroids based on the mean po- 
sition of member galaxies or X-ray flux have off- 
sets from other center definitions that are roughly 
consistent with their larger statistical uncertainties 
(~ 50- 150 kpc). 

2. The offsets between centers produce a measurable 
signal in the lensing profile. Stacking the signal 
around a bright or massive galaxy tends to produce 
a larger signal at small scales than any centroid 
not located on a galaxy, and the difference in these 
signals is greater than expected from the stellar 
mass of the galaxy. 

3. Among the candidate centers we have tested, the 



brightest or the most massive galaxy near the X-ray 
centroid appear to be the best tracers of the center 
of mass of halos. Centering definitions based on 
the centroid of member galaxies have larger offsets 
and uncertainties. 

4. Groups that have ambiguous centers because of 
multiple bright or massive galaxies have lensing 
signals that suggest a lower mass than expected 
given their X-ray luminosity and in some cases ap- 
pear disturbed. These are possibly merging sys- 
tems, and the property of having discrepant candi- 
date centers gives a simple observational indicator 
to identify them. 

5. In groups with a clear central galaxy, offsets be- 
tween the galaxy and the halo center are fairly 
small (< 75 kpc). The offset is somewhat degener- 
ate with the amount of substructure around the 
galaxy and with the concentration of the group 
halo. 

These findings apply to our group sample but our ap- 
proach can readily be applied to other group and cluster 
data sets, and to different analyses such as satellite dy- 
namics and richness-based mass estimators. Given the 
level of disagreement among our candidate centers, we 
advise testing different centers to determine the degree to 
which centering choices affect a given analysis. Our find- 
ing that groups with ambiguous centers are less massive 
for their X-ray luminosity or have disturbed mass pro- 
files suggests that these groups could be identified and 
excluded if an analysis calls for a clean sample of halos. 
Additionally, a probabilistic approach to centering algo- 
rithms could provide information about the confidence 
in a given center allowing for appropriate weighting of 
different systems. 

Much larger samples of groups and clusters are being 
constructed with ongoing and upcoming surveys such as 
the South Pole Telescope and Atacama Cosmology Tele- 
scope, the Dark Energy Survey, and eROSITA. This ini- 
tial study with a modestly-sized sample of groups bene- 
fits from high-resolution X-ray selection which provides a 
good starting point for finding centers, but a similar ap- 
proach can be applied to optical or SZ-selected catalogs. 
With larger samples, we can hope to improve constraints 
on the density profiles of dark matter halos including 
their concentration and inner slope, as well as the mass 
distribution in subhalos and the effects of mergers. Ad- 
ditional constraints on the mass distribution from strong 
lensing or stellar kinematics can provide interesting con- 
straints on a range of scales and improve models for the 
distribution of offsets between galaxies and halo centers. 
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